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I. INTRODUCTION 



(Dated: March 11, 2008) 

We formulate a bosonic dynamical mean-field theory (B-DMFT) which provides a comprehen- 
sive, thermodynamically consistent framework for the theoretical investigation of correlated lattice 
bosons. The B-DMFT is applicable for arbitrary values of the coupling parameters and tempera- 
ture and becomes exact in the limit of high spatial dimensions d or coordination number Z of the 
5_l ■ lattice. In contrast to its fermionic counterpart the construction of the B-DMFT requires different 

scalings of the hopping amplitudes with Z depending on whether the bosons are in their normal 
state or in the Bose-Einstein condensate. A detailed discussion of how this conceptual problem 
can be overcome by performing the scaling in the action rather than in the Hamiltonian itself is 
presented. The B-DMFT treats normal and condensed bosons on equal footing and thus includes 
the effects caused by their dynamic coupling. It reproduces all previously investigated limits in 
parameter space such as the Beliaev-Popov and Hartree-Fock-Bogoliubov approximations and gen- 
eralizes the existing mean-field theories of interacting bosons. The self-consistency equations of the 
B-DMFT are those of a bosonic single-impurity coupled to two reservoirs corresponding to bosons in 
the condensate and in the normal state, respectively. We employ the B-DMFT to solve a model of 
itinerant and localized, interacting bosons analytically. The local correlations are found to enhance 
the condensate density and the Bose-Einstein condensate (BEC) transition temperature Tbec- This 
, effect may be used experimentally to increase Tbec of bosonic atoms in optical lattices. 

PACS numbers: 71.10.Fd, 67.85.Hj 

£ 
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The observation of Bose-Einstein condensation (BEC) in ultra-cold, atomic gases has greatly stimulated research 
into the properties of this fascinating quantum state of matter i In particular, experiments with alkali atoms confined 
I> in optical lattices 2 -^ have renewed the theoretical interesl^LS in the physics of strongly correlated bosons on lattices, 
ON , which promises significant new insights and even applications in fields such as quantum computing^ 

The investigation of correlated lattice bosons is not only relevant for ultra-cold bosonic atoms in optical lattices, 
but has a long history starting with Matsubara and Matsuda' o 10 ' n i 12 formulation of a lattice model of liquid 4 He. 
A new direction of research was initiated by Fisher et a/.,™ who studied lattice bosons with and without disorder 
to explore the superfluid-insulator transitio n 14 ' 15 i 16 i 17 i 18 and boson localization observed in 4 He absorbed in porous 
media. 19 Granular superconductors forming weak Josephson junctions have also been described by interacting lattice 
bosons. 20 ' 21 Recently quantum phase transitions in magnetic systems such as TlCuClaj^ 2 . which can be induced by 
tuning the magnetic field, have been interpreted as the BEC of magnons. 2 ^. Bosonic supersolid a 24 ' 25 promise to be yet 
another fascinating state of bosonic matter. 

In this paper we formulate the first comprehensive, thermodynamically consistent theory of correlated lattice 
boson systems, namely a bosonic dynamical mean-field theory (B-DMFT) which is applicable for arbitrary values 
of the coupling parameters and temperature. The B-DMFT includes all local, dynamical correlations of the many- 
boson system and becomes exact in the limit of infinite space dimensions in analogy with its successful fermionic 
counterpart i 26 ' 27 i 28 ' 29 With the B-DMFT we are able to solve a lattice model of itinerant and localized, interacting 
bosons. The local correlations lead to an enhancement both of the BEC transition temperature Tbec and the 
condensate fraction as compared to the non-interacting system. Hence bosonic correlations can be employed in the 
laboratory to reach higher values of Tbec- 

This paper is organized as follows: In Sec. II we introduce the bosonic Hubbard model and explain the spe- 
cific problems arising in the construction of a bosonic dynamic mean-field theory (B-DMFT) in the limit of large 
coordination number Z of the lattice, namely, the problem of how to scale the hopping amplitudes with Z. The 
self-consistency equations and the general structure of the B-DMFT are discussed in Sec. III. The comprehensive 
nature of the B-DMFT is demonstrated in Sec. IV by explicitly reproducing results previously obtained in special 
limits of parameter space and by deriving other bosonic mean-field theories. In Sec. V the B-DMFT is employed to 
solve a bosonic version of the Falicov-Kimball model and it is shown that correlation effects lead to an enhancement 
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of Tbec- Conclusions and an outlook in Sec. VI close the presentation. 



II. CORRELATED LATTICE BOSONS 



A. 



Generalized bosonic Hubbard model 



In the following we consider a many-particle system with different species of bosons as it can be realized in optical 
lattices^! 4 - This may either involve different atoms as, for example, in a binary mixture of 87 Rb and 7 Li, or one type 
of atom in different hyper-fine quantum states such as 87 Rb, where the total nuclear spin 7 = 3/2 adds to the spin 
S = 1/2 of the valence s-electron, giving states with f = / + 5 = 1 or 2. Such systems may be modelled by a 
generalized bosonic Hubbard Hamiltonia n 5 ' 10 ' 13 ' 30 



where 7ij„ = b iv bi V is the occupation number operator for bosons of species v. Furthermore, are hopping amplitudes 
of i^-bosons and are local density-density interactions between \i- and ^-bosons on the same lattice site. An 
exchange interaction for spinor bosons can be easily included. In general the many-boson model ([T]) in unsolvable. 



The explanation of experiments with correlated lattice bosons in quantum optics and condensed matter physics 
requires a comprehensive theoretical scheme for the investigation of the Hamiltonian (fTJ). In particular, it must be 
capable of describing thermal and quantum phase transitions and thus provide the phase diagram and the thermody- 
namics for the entire range of microscopic parameters. In the case of lattice fermions such a framework already exists: 
the dynamical mean-field theory (DMFT)2&. Indeed, the DMFT has proved to be a very successful, comprehensive 
mean-field theory for models and materials with strongly correlated electrons! 29 ' 31 In particular, it provides a quanti- 
tative description of the Mott-Hubbard metal-insulator transition, photoemission spectra, magnetic phases, and other 
correlation induced phenomena. The DMFT has the virtue of becoming exact in the limit of infinite space dimensions 
d or, equivalently, infinite coordination number Z, i.e., number of nearest neighbors (Z — Id for a d-dimensional 
hypercubic lattice)^ This limit is well-known to produce mean-field theories which are diagrammatically controlled 
and whose free energies have no unphysical singularities (e.g., the Weiss mean- field theory for the Ising or Heisenberg 
spin models). 32 To obtain a physically meaningful mean- field theory the free energy of the model has to remain finite 
in the limit d or Z — > oo. 2 - This requires a suitable scaling of the coupling parameters with d or Z, e.g., J — > J /Z, J— 
const., for Ising spins with nearest-neighbor coupling J. While for the Ising model the scaling is self-evident this is 
not so for more complicated models. Namely, fermionic or bosonic many-particle systems are usually described by a 
Hamiltonian consisting of several non-commuting terms each of which is associated with a coupling parameter, e.g., a 
hopping amplitude or interaction. In such a case the question of how to scale these parameters has no unique answer 
since this depends on the physical effects one wishes to explore! 33 ' 34 In any case, the scaling should be performed such 
that the model remains non-trivial and its free energy stays finite in the Z — > oo limit. By "non-trivial" we mean 
that not only (Ho) and (iJj nt ), but also the competition between these terms as expressed by ([Hq, _ff; n t]), should 
remain finite; here (...) denotes the quantum and statistical average of operators. In the literature on lattice bosons 
the d — ► oo limit was so far considered only in connection with the distance- independent ("infinite-range") hopping 
of the boson s 15 ' 17 ' 18 in which the mean-field theory of Fisher et al^ for the Bose-condensed phase becomes exact. 
As will be discussed below this is a static mean-field theory since normal and condensed bosons are not dynamically 
coupled. In particular, in the normal phase one has ([-ffn, -Hint]) = and the lattice problem is reduced to a single-site 
("atomic") problem where particles are immobile. Another static mean-field theory is the Bogoliubov approximation, 
which yields a good weak-coupling mean-field theory for bosons in a continuum. For lattice bosons this approximation 
fails to describe the Mott superfluid-insulator transition. 

Evidently, a dynamical mean-field theory is needed to describe the rich physics of interacting lattice bosons, e.g. 
cold atoms in optical lattices? 2 - within one conceptual framework. A comprehensive DMFT for correlated lattice 
bosons, i.e., a theory which can describe normal and condensed bosons on the same footing, did not exist up to now. 
In the following we will discuss the conceptual problems which prevented the formulation of such a bosonic DMFT, 
and how they can be overcome. 




(1) 



B. Construction of a comprehensive mean- field theory 
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C. Lattice bosons in infinite dimensions: Different scaling for Bose-Einstein condensed and normal bosons 



A macroscopically large number of bosons can condense into a single quantum state. This BEC may be detected 
in the spectral decomposition of the one-particle density matrix 

Pij = ( b l b j) = E A "Oi"> ( 2 ) 
a 

where b\ and bi are creation and annihilation operators, respectively, for a boson at a lattice site i, with <f>i a as the 
corresponding wave function. For simplicity we discuss here only a single species of bosons so the index v can be 
omitted. When BEC occurs one of the eigenvalues becomes macroscopically large, Ao = No ~ O(N), where TV is the 
total number of bosons. The density matrix then decomposes into 

Pij = N </>* <f> j0 + Pij, (3) 

where the second term corresponds to non-condensed, "normal" bosons. The first term has the remarkable property 
that it does not decrease even at large distance — ||Rj — Rj|| between the bosons at sites i and j. Here ||R|| 
denotes the length of R obtained by counting the minimal number of links between two sites on a lattice. By contrast, 
the second term in p^ decreases with increasing R y . 

This has immediate consequences for the kinetic part of the Hamiltonian 

H = -tJ2 b\b 3 (4) 

<ij> 

with —t as the amplitude for hopping between nearest neighbor sites i and j. For a uniform BEC with density 
no = N /Nl, where Nl is the number of lattice sites, one has </>* o 0jo = l/^Vt such that the kinetic energy is given by 



Skin - ~t E Pv ~ ^kin + ^£n rmal ! ( 5 ) 

where 



<ij> 



s k B in c = -t E ^o, s£°r al = -*E^'- ( 6 ) 

<ij> <ij> 

To derive a mean-field theory for lattice bosons via the limit of high spatial dimensions the energy density E^/Nl 
needs to remain finite for d or Z — > oo. Since the energy density of the condensate, E^ c /Nl — Ztno, is proportional 
to Z a non-trivial limit Z — > oo is obtained only if the hopping amplitude is scaled as t = t/Z, with t — consti 17 i 18 In 
the case of normal lattice bosons (or fermions) the situation is characteristically different. Since pij is the transition 
amplitude for the hopping of a boson from j to one of the Z neighboring sites i the respective normalized hopping 
probability is \pij\ 2 cx 1/Z, whence pij ~ \ j\[Z. For the energy density of the normal bosons, E£°* mal /Nl oc Ztptj, 
to remain finite for Z — > oo the hopping amplitude must therefore be scaled as in the case for fermions, namely, as 
t = t/ v^.— In the more general case of hopping between sites i and j which are not nearest neighbors, the amplitudes 
tij have to be scaled as 

Uj = tij/(z R »y, (7) 

where s = 1 ("integer scaling") if the bosons are quantum condensed and s = 1/2 ("fractional scaling") if they are in 
the normal state. 

The total energy of a single species (v = 1) of correlated lattice bosons described by the Hamiltonian ([1]) is given 

by 

E = - f E N °-t E + \ u zZ^ ni - !))• ( 8 ) 

<i,j> <i,j> i 

If the scaling of the hopping amplitudes in the limit Z — > oo is performed on the level of the Hamiltonian (or the 
energy E) two cases have to be distinguished: 

(i) N = 0: In the absence of a BEC fractional scaling (t — t/yZ) has to be employed to arrive at a finite value 
of E for Z — > oo. We note that the interaction is purely local and hence independent of the spatial dimension of the 
system; consequently U need not be scaled at all. 
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(ii) A*o 7^ 0: In this case integer scaling (i = t/Z) has to be employed. Thereby the contribution of the condensate 
(first term in ([8]) ) remains. However, at the same time the contribution of the non-condensed (normal) bosons to the 
kinetic energy is suppressed oc \j\[Z. The normal bosons thus become immobile. It can be shown rigorously that 
the mean-field equations obtained in this way are equivalent to those derived by Fisher et alr^ which are known to 
be exact in the limit of infinite-range hopping (see Section IV. C). 

The above discussion shows that the construction of a mean-field theory via the limit Z — > oo which is based on a 
scaling of the hopping amplitudes in the Hamiltonian is either restricted to the normal state, or removes the normal 
bosons from the problem. This is unsatisfactory since the important dynamical coupling of normal and condensed 
bosons is then eliminated from the outset. During the last 15 years the problem of how to scale the hopping amplitudes 
without eliminating the normal bosons presented an unsurmountable obstacle for the formulation of a bosonic DMFT. 
Indeed, our discussion shows that such a theory cannot be formulated on the level of a fully microscopic Hamilton 
operator, i.e., without making an additional Bogoliubov mean- field type assumption. 35 

At this point it should be pointed out that there exists no a priori condition according to which the scaling of the 
hopping amplitudes has to be performed in the Hamiltonian. Indeed, for the free energy of the model to remain finite 
in the limit Z — > oo the scaling can equally be performed in the partition function (or in the action entering in the 
functional integral which determines the partition function), from which the free energy is calculated. 

III. BOSONIC DYNAMICAL MEAN-FIELD THEORY (B-DMFT) 
A. General structure of the B-DMFT 

We will now show that the long-standing problem of the scaling of the hopping amplitudes with the coordination 
number Z or the dimension d can be resolved by considering the large dimensional limit not in the Hamiltonian but in 
the action determining the Lagrangian density. Namely, integer scaling is employed whenever a hopping amplitude is 
associated with anomalous expectation values (6,-) and (b\), while fractional scaling is used otherwise. The B-DMFT 
obtained in this way treats normal and condensed bosons on equal footing and is thus able to describe both phases, 
including the transition between them, in a thermodynamically consistent way. 

The general structure of the B-DMFT is shown in Fig. ([1]). In the limit d — > oo the bosonic many-body lattice 
problem is mapped onto a single-site problem with integer occupation. This site is coupled to two particle reservoirs, 
one representing normal the other quantum condensed bosons. These reservoirs represent the Weiss-type molecular 




time 



FIG. 1: Bosonic dynamical mean-field theory (B-DMFT): Within the B-DMFT the full many-body lattice problem is reduced 
to a single-site problem which is coupled to two reservoirs corresponding to bosons in the Bose-Einstein condensate and in 
the normal state. The integer occupation of the site changes in time and is determined by the local interactions and the 
time-dependent properties of the particle reservoirs. Although the total number of bosons is preserved, particles are scattered 
between the normal and the condensate reservoirs via the single site as shown by arrows. This schematic picture visualizes the 
idea of DMFT for lattice bosons in analogy to the fermionic counterpart described in Ref. [29l . 
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fields of the B-DMFT. Their properties are determined self-consistently by the B-DMFT equations. Particles hop 
onto and off the site as a function of time, thus changing the total number of bosons of the reservoirs. Therefore local 
correlations lead to a dynamical depletion or filling of the condensate. 

B. B-DMFT equations 

Here we present the self-consistency equations of the B-DMFT for the model (fTJ) ; their derivation is discussed in 
detail in Appendix A. The time evolution of i/-bosons on the single site i = is represented by the local propagator 
(Green function) 

G v (T) = -(T T K(T)bt(0))s^, (9) 

where we used the imaginary time, finite temperature formalism and Nambu notation bj, = (&*,&„); T T is the time 
ordering operator. The diagonal elements of the Green function matrix in Nambu space represent the quantum- 
mechanical probability amplitude for creating a boson on a site at one particular time and annihilating it after a time 
t, or the time inverted process. The off-diagonal elements, present only in the BEC phase, represent the amplitudes 
for creating or annihilating two bosons at different times. In the path integral formalism the probabilities of such 
processes are determined by the local action, which in the B-DMFT is given by 

Sioc = - f dr f dT'Y^hi{T) g-\ T -r')K(r')+ fdrV^(T)K(r)-g+ f dr V) kM (r)b v (r](10) 
J o Jo u Jo M „ 1 Jo 

Here k v is a numerical factor depending on the lattice structure, i.e., k v = J2i^o ^io/^ Rio f° r d — > oo, and k v = J^i^o tfo 
for an approximation in finite dimensions. The free ("Weiss"—) mean field propagator Q~ , which is determined by 
the properties of the reservoir of normal i/-bosons, is related to the local propagator G„ by the Dyson equation 

Gv l {iu n ) = G" 1 + £„(iu; n ) = (w n <Ta - n„l) - A v (iui n ), (11) 

where Q, and X arc also matrices in Nambu space. Here u> n = 2im/ (3 are even Matsubara frequencies with the inverse 
temperature = 1/ksT and Ti^{iuj n ) is the momentum-independent (local) self-energy. The quantity A„ describes 
the resonant broadening of quantum-mechanical states on a lattice site and may be interpreted as a hybridization 
of bosons on that site with the surrounding bosonic bath. This hybridization function is determined by the local 
correlations through eq. (fT0|) . The third term in (fTTJ|) describes the coupling of a local boson to the condensate, the 
latter being represented by an order parameter 3>J,(t). In our formulation this term arises naturally in the case of 
BEC and does not require a Bogoliubov substitution^ 

The second B-DMFT equation is given by the lattice Hilbert transform 

G v (iu} n ) = 22 [(iw n <T3 - (e£ + - S l /(«w„)]" 1 , (12) 

k 

where e£ is the dispersion relation of non-interacting z/-bosons, \x v the chemical potential, k the wave vector, 1 a unity 
matrix, and 03 the Pauli matrix with ±1 on the diagonaL 27 ' 28 Eqs. (|9I12[) are the counterparts to the self-consistency 
equations of the DMFT for correlated lattice fermions. However, here these equations contain the condensate wave 
function i.e., the order parameter of the BEC, which enters as a source field in the action ©. It can be determined 
exactly by calculating the average 

*„(r) = (b u (r)) Sloc (13) 

together with eqs. (J9j) and (Tl2|) . We note that in equilibrium the time dependence enters via a trivial exponential 
factor e~^ vT which can be eliminated by a gauge transformation^ 

Eqs. (|9l fl3| constitute the B-DMFT solution of the generalized bosonic Hubbard model (fTJ). These equations are 
exact in the d — > 00 limit and provide a comprehensive, thermodynamically consistent and conserving approximation in 
finite dimensions. In other words, the B-DMFT derived here is the first mean- field theory for correlated lattice bosons 
which has all the attractive features characterizing the now well-established fermionic DMF T 28 ' 29 In particular, the 
B-DMFT can be expected to be the best approximation to many-boson problems with strong local correlations since 
the on-site quantum fluctuations of the spin or density are treated exactly. Spatial correlations are neglected but can 
be restored, e.g., within cluster extensions of the B-DMFT. Furthermore, long-range ordered phases can be described 
within the B-DMFT by properly choosing the self-consistency conditions in analogy with the fermionic case^ 
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C. B-DMFT and Gross-Pitaevskii equation 

The exact Euler-Lagrange equation of motion for the field b(r) is obtained from the stationary conditions 
SS\ oc [ bl, b v ]/S bl = of the local B-DMFT action and is given by 

dr K{t) - f dr'[ Al\r - r') b v (r') + A* 2 (r - r') &t( T ')] + k v b v (r) + £ &t (r) b^r) b u (r) = \i v b u (r).(U) 
J ° „ 

If one replaces each field b(r) by its expectation value (the order parameter <& v ) one arrives at 

d T $„(r) - / dr'l Delta? (t - r') + A^ 2 (r - r') $,>')] + «„$„(,-) ^(t)| 2 *„(t) = ^ ^(rX15) 

This is a generalization of the standard Gross-Pitaevskii mean-field equation, a non-linear differential equation for the 
condensed bosons which can be derived within the time-dependent Hartree-Fock approximation. 35 At present it is not 
clear whether the replacement b v (r) — > $ u (t), i.e., the factorization of the correlation function ( 6t,(r) 6 M (r) b u {r)) = 
| "^(t)! 2 $j/(t), holds rigorously as in other mean-field theories i 47 ' 48 ^ 49 ' 50 Eq. (fl~5|) is the classical equation of motion 
of the condensate of lattice bosons in the d — > oo limit. The second term on the l.h.s. of (fl~5|) describes the retarding 
effects of normal bosons on the condensate and makes the generalized Gross-Pitaevskii equation a non-linear integro- 
differential equation. We note that in the standard Gross-Pitaevskii equation the hybridization term A„ is missing. 
Once A y has been determined by solving the B-DMFT self-consistency equations ([T5|) can be used to determine 
3>„(t) and then calculate any response function of the condensate due to external perturbations, e.g., to describe 
Bragg scattering. In equilibrium & v can be expected to be independent of r. Then the stationary solution of (|15p is 
easily obtained by solving a set of linear equations. For example, for the spinless bosonic Hubbard model (no index 
v) we obtain |$| 2 = [fi — n + A 11 (w„ = 0) + A 12 (w„ = 0)]/U. The condensate density depends explicitly on the zero 
mode components of the hybridization functions for the normal subsystem. 

IV. THE B-DMFT IN DIFFERENT LIMITS OF PARAMETER SPACE 

The B-DMFT is a comprehensive mean-field theory for correlated bosons on a lattice, meaning that the theory is 
valid for all input parameters and temperatures. This is an essential prerequisite for obtaining a reliable, approximate 
description of those parts of the phase diagram which cannot be studied perturbatively. Above all, the B-DMFT 
reproduces all known results obtained in special limits of the parameter space as depicted in Fig. [2] In particular, the 
exactly solvable limits of free and immobile bosons, respectively, and of well-known static mean-field approximations 
can be obtained directly from the B-DMFT. In the following we discuss these limits in details. 

A. Free bosons 

In the non-interacting limit, = 0, the problem is trivially solvable in all dimensions. Since there is no interaction, 
all correlation functions factorize and the cavity method employed in Appendix A becomes exact. In this case the BEC 
is described within the grand-canonical ensemble by a non- vanishing order parameter <E>. However, the off-diagonal 
Green (hybridization) functions of the normal bosons are zero. Explicitly, the local action has a bilinear (Gaussian) 
form 

^on-interacting = _ t* ^ /"^ ^b^r) ^( T -/) ^ + dr £ (r)b „ ( r)) (1 6 ) 

Jo Jo Jo 

with Q~ x (t — t') = (— d T <T3 + fj, v l)6(T — t') — A„(t — r'). All known equations, e.g., that for the particle number 
or the compressibility, can be easily derived from this action (see Appendix B). We note that in the non- interacting 
limit the hopping amplitudes need not be scaled for d < oo. However, to obtain a meaningful limit d — ► oo the scaling 
scheme introduced in this paper is necessary. Otherwise the condensate and normal bosons would not be treated on 
equal footing. That is, if only fractional scaling is employed, one obtains spurious infinities in the condensate phase, 
whereas if only integer scaling is used, the normal bosons (which contribute significantly at temperatures close or 
above Tbec) become immobile. 
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Bosonic Hubbard Model (t»,U) 
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s = \ for condensed bosons 
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FIG. 2: Relation of the B-DMFT to other approximations and exact limits; see Sec. IV. 



B. Immobile bosons ("atomic limit") 

In the atomic limit, t£j- = 0, all lattice sites are decoupled and the particles are immobile. In this case the order 
parameter = since no condensation is possible, and also the hybridization function A„ = 0. For arbitrary 
dimensions the exact action is then given by a sum over all equivalent sites with the same local action 



Si] 



(17) 



: = - f dT Y, h l( T ) + 1*1) Mr) f(iTr^(f)K(T)-y, 

as obtained from the action (fT0|) with <& v = and A„(r) = in Qy X (r — t') = (— <9 T <73 + h v 1)6(t — t') — A„(r — t'). 

C. Mean-field theory of Fisher et al. 

The mean-field field theory of Fisher et al*£ is known to be the exact solution of the bosonic Hubbard model when 
the hopping amplitude is independent of distance and is scaled with the number of lattice sites Nl, i.e. = i/N n 13 i 41 
This is also called the limit of "infinite-range hopping" , 13 i 41 i 42 The free energy density of the bosonic Hubbard model 
in this limit has the form 



= F at - kTln(T T e- f« dT ^^ b »^) Hat - V^|^| 2 . 



(18) 



Here the atomic part F at — — fcTlnTr exp(— f3H at ) is given by the Hamiltonian H at , which is obtained from the lattice 
Hamiltonian by setting all hopping amplitudes equal to zero. The average in the second term is taken with respect to 
H at . The stationarity condition for the free energy (fT5|) with respect to i.e. QF ln ^ mt °-^& c j Q^ u = rj then yields 
the self-consistent mean-field equation 



,*„ = kT 



(T T e-/. s *E»^*lW{r) )fftt 



(19) 



Since the normal bosons are immobile, and thus not dynamically coupled to the condensed bosons, the theory of 
Fisher et al.^ is a static mean-field theory. 
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We now show that the mean-field equation ([19]) can also be obtained in the d — > oo limit as put forward in 
Refs. [l5llT7lfl8l . Indeed, by employing the cavity method for lattice bosons and using only integer scaling for the 
hopping amplitudes the local action in the d — ► 00 [Z — > 00) limit takes the form 

s£r°' intcgcr scalins = - / dT E b -( T ) + b -( r ) 

+ / dr^%^(r)K(r)-^]+ / dr5>„*t ( r ) bl/ (r). (20) 

This expression differs from the local action in the atomic limit, ()17|) . by the presence of the last term which describes 
the condensate. However, in equilibrium the BEC order parameter is time independent in which case (|20p yields the 
free energy density as 

pdr+oo, integer scaling = ^ _ fcy^ exp ^_ ^ dr £ ) H „ . (21) 

The BEC order parameter & v obeys the self-consistent equation 

*„ = —. 4 r- ' / D & b *l b - e ~ s;ir "' '° teSM Bcalta8 [ fe ' b * ; *"]. (22) 

r7fl — >oc, integer scaling / L > J v / 



loc 

Other correlation functions and observables can be determined similarly. 

The free energy density (f2"Tj) is seen to differ from (fTS)) only by the absence of the last term proportional to the 
density of the condensate. Nevertheless, the equations for the BEC order parameter as well as correlation functions 
and observables are the same. Thus the approximation of constant hopping amplitude ("infinite-range hopping") and 
the d — ► oo limit with integer scaling give rise to the same mean-field equations. At T = these equations can also be 
derived by yet another approximation, namely a variational method using a Gutzwiller-type wave function. 15 - 43 The 
mean- field theory of Fisher et al. 13 and its generalization to spinful bosons were widely used to investigate quantum 
phase transitions and the phase diagrams of correlated lattice boson systems and of mixtures of lattice bosons and 
fermions^ 5 -^^ 5 .^ 

Eq. (|21[) and l|22p can also be obtained directly from the B-DMFT self-consistency equations by neglecting the 
hybridization function, i.e. by setting A„ = 0. Then the local action of the B-DMFT, (fTT))) . is the same as that in 
([20]) . Furthermore, (fT2|) is satisfied automatically in this limit since only the state with k = is taken into account. 
It should be noted, however, that in the absence of the hybridization function A„ the non-interacting limit of the 
normal bosons cannot be reproduced, i.e., the mean-field theory of Fisher et ali^ does not describe the limit of free, 
normal bosons. 



D. Weak-coupling (Bogoliubov) mean-field theory 



A perturbation expansion to first order in is equivalent to the Hartree-Fock-Bogoliubov approximation with the 
static self-energy E 11 = 2 ^ U uu n^ G - 2 U^G„(uj n )/p and S^ 2 = £ M U^^ EC . For such a self-energy the 

self-consistency condition (| 1 2[) is equivalent to the self-consistent Hartree-Fock-Bogoliubov approximation (sometimes 
called "first-order Popov" approximation) ! 35 i 36 This self-consistent approximation is known to lead to a gapped 
spectrum in the condensed phase because off-diagonal elements in the self-energy are calculated in higher order due to 
self-consistency. By contrast, the Bogoliubov approximatio n 35 ! 36 is obtained if only particular diagrams corresponding 
to the self-energies E 11 = 2^^ t/^n™ and S 12 = ^ U^^n^ EC are taken into account. 

The second-order expansion contains many diagrams, see Refs. |35ll36l . Checking term by term we find that the 
B-DMFT reproduces the Beliaev-Popov approximation (sometimes called "second-order Popov approximation" ^ if, 
in addition, in the latter approximation only local irreducible self-energy diagrams (consistent with the d — > oo limit) 
are retained. 



V. B-DMFT SOLUTION OF THE BOSONIC FALICOV-KIMBALL MODEL 



We now apply the B-DMFT to study BEC in a mixture of two different species of bosons: itinerant 6-bosons and 
immobile /-bosons. We assume the 6-bosons not to interact with each other but only with /-bosons, while f-bosons 
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interact also mutually, i.e., Ubb — 0, Ubf > 0, lift > in JT]). We call this the bosonic Falicov-Kimball (BFK) 
model since it is a bosonic generalization of the Falicov-Kimball model for fermions which has been widely studied in 
condensed-matter physics^ 7 . Experimentally such a system can be realized by loading an optical lattice either with a 
mixture of two different species of bosonic alkali atoms (e.g., 7 Li and 87 Rb), or by one kind of atom with two different 
hyper-spin states (i.e., with F = 1,2 and specific values of the z-component of F). In addition, the electric fields 
generating the potentials of the optical lattice and the external magnetic field controlling the Feshbach resonances 
should be tuned such that one species of particles is immobile and the other is non- (or only weakly) interacting. 
The realization of a fermionic Falicov-Kimball model by cold fermionic atoms in an optical lattice was discussed in 
Refs. HM 

It is important to note that, in spite of the immobility of the /-bosons, the BFK model is still a many-body 
problem because the immobile particles are thermodynamically coupled to the mobile particles by the interaction. 
In particular, the optimal configuration of the localized bosons depends on the interaction, temperature and density 
of the particles. In the fermionic counterpart one finds that the position of the immobile particles is either random 
or long-range ordered; phase separation between these two components can also occur. The numerical solution 
of the Falicov-Kimball model is limited to small lattices and requires an annealed average over a large number of 
configurations of immobile particles. 

For the BFK model the local impurity problem can be integrated analytically. The self-consistency equations can 
be then solved by standard numerical techniques. Since the /-bosons are immobile, their number on each site is 
conserved. Hence, the /-boson subsystem cannot undergo BEC and the occupation number operator nf of the single 
site becomes a classical variable with rif = 0, 1, 2, ... The local action © is then quadratic in the bosonic operators. 
Consequently, the local propagator Gb{iu n ) for 6-bosons and the local partition function Z\ oc (nb, fif), and thereby 
the BEC transition temperature Tbec for the 6-bosons, can be evaluated directly. The local partition function of the 
BFK model is determined by 

oo v n 

Zioc(M6,i«/)= £ e^^-^^ZU^-UbfUf,^), (23) 

n/ =0,1,2,.. 

where 



^° c (M 6) M/)~e ^WJJU— ^7-^ (24) 

is the partition function for Ubf = 0. The local propagator for normal 6-bosons is given by 

oc 

G b (iu> n ) = ~— JT 1 aTT^" ( 25 ) 

iu n + fib - UbfUf - A b {iuj n ) 



n/=0,l,2,. 

Here 

Zi oc (fi b , fif) 

is the probability for the single site to be occupied by exactly rif = 0,1,2,... bosons. For hard-core /-bosons 
(Uff = oo) this leads to w nf= \ = hf and w nf= a = 1 — hf. In this case hf rather than ixj is used as an independent 
thermodynamical variable. The propagator (|25[) describes quantum and thermal fluctuations of normal bosons. In 
the absence of the interaction between 6-bosons, off-diagonal terms in the local propagator are zero. The Gross- 
Pitaevskii equation is then obviously exact and reduces to a homogeneous, linear equation of the form [— iuj n — Lib + 
Kb + Ab(iuj n )](j>b('i^n) = for each Fourier component. For oj n ^ the only solution is </>&(cJ n 7^ 0) = 0. The static 
{uj n = 0) component of the BEC order parameter is finite if fib = k& + A&(0) and must be determined by fixing the 
average density of 6-bosons. 

A striking result obtained for this model is an enhancement of Tbec f° r increasing repulsion between the 6- and 
/-bosons, with a maximum of Tbec a t intermediate values of Ubf and a saturation at large Ubf- This behavior is 
explicitly seen in Fig. [SJ'V for hard-core /-bosons {Uff = 00, i.e., nf — 0, 1) on a simple cubic lattice where we plotted 
the relative change of Tbec with respect to Tg EC in the non-interacting system. The increase of Tbec is due to the 
blocking of a fraction of sites by heavy atoms which increases the density of the 6-bosons. However, this argument 
cannot explain the non-monotonicity of Tbec vs Ubf shown in Fig. [3^. In fact, the maximum is due to the correlation 
induced band splitting and the narrowing of the lower subband (see Fig. [3p) which lead to an increase and decrease 
of TbeCi respectively. Furthermore, at fixed temperature T the average condensate density nf EG (T) is also found to 
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FIG. 3: Bose-Einstein condensation of a mixture of itinerant and localized, correlated lattice bosons: (A) Enhancement of the 
BEC transition temperature with increasing interaction strength Ubf in a two component boson mixture with different densities 
n/ of the localized /-bosons. (B) Dependence of the condensate fraction nf EC (T)/n{, on temperature for different interactions 
Ubf at fif — 0.8. (C) Spectral functions for different values of Ubf at fif = 0.8. The increase of Tbec and the condensate 
fraction with increasing Ubf and fif is caused by correlation effects leading to a redistribution of the spectral weight for the 
6-bosonic subsystem. The correlation gap opens when Ubf exceeds a critical value which depends on fif. The opening of the 
gap is not associated with a phase transition of the mobile bosons. Results are obtained for a three-dimensional cubic lattice 
with unit band-width and fib = 0.65. In the hard-core limit the spectral functions are temperature independent because the 
occupation probability of /-bosons is either fif or 1 — fif. 



increase with increasing repulsion (see Fig. 03), although the interaction induced scattering between bosons usually 
removes particles from the condensate, thereby reducing its density^ At zero temperature all &-bosons are in the 
condensate. Similar results are found for other lattices, Cf. Fig. 0] 

These results originate from local correlations which are captured exactly by the B-DMFT, but not by conventional 
approximations. The consequences can be inferred by considering the total density of 6-bosons 

n b = n b (T) + / du — — -. (27) 

J exp(cj/T) - 1 

The second term gives the contribution of normal &-bosons for which the spectral function A b (uj) = — ImG b (uj)/ir is 
shown in Fig. [3p for different U b f values. The spectral weight is seen to be strongly redistributed, forming lower- 
and upper-Hubbard subbands at low and high energies u, respectively, which are separated by the energy Ubf- We 
note that the splitting and rounding of the shapes are genuine correlation effects. The occupation probability of 
normal 6-bosons, i.e., the Bose-Einstein distribution function, decays exponentially with increasing lo. This implies 
an extremely small particle number in the upper Hubbard subband at large U b f. Since the total number of 6-bosons 
is constant, the particles are necessarily transferred into the condensate. Hence, at fixed temperature fif increases, 
thereby enhancing Tbec- The spectral weight contributing to the upper Hubbard subband is also proportional to fif, 
implying that Tbec increases with fif, too. 

The B-DMFT prediction of an increasing Tbec an d condensate density due to local correlations are expected to 
be observable in mixtures of mobile and localized bosons on three dimensional optical lattices. Such correlations can 
thus be employed in the laboratory to enhance Tbec of bosonic condensates. We also note that on bipartite lattices 
with special densities of bosons, e.g. fif = n b — 0.5, long-range order in the density of the /-subsystem and, in turn, 
a supersolid phase in the 6-subsystem is expected to form. Obviously the physics of this seemingly simple bosonic 
model is extraordinarily rich. 
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FIG. 4: Bose-Einstein condensation of a mixture of itinerant and localized, correlated lattice bosons: In contrast to Fig. 2 
results here are obtained for a Bethe lattice with infinite Z, unit band-width, and fib = 0.5. They are exact for the bosonic- 
Falicov-Kimball model. (A) Enhancement of the BEC transition temperature with increasing interaction strength Ubf in a two 
component boson mixture with different densities n/ of the localized /-bosons. (B) Dependence of the condensate fraction 
nf EC (T) /fib on temperature for different interactions Ubf at ftf = 0.8. (C) Spectral functions for different values of Ubf at 
fif = 0.8. The increase of Tbec and the condensate fraction with increasing Ubf and n/ is caused by correlation effects leading 
to a redistribution of the spectral weight for the 6-bosonic subsystem. 



VI. CONCLUSIONS 



In this paper we derived the first comprehensive, thermo dynamically consistent theoretical framework for the 
investigation of correlated lattice bosons — a bosonic dynamical mean- field theory (the B-DMFT). In analogy to its 
fermionic counterpart the B-DMFT becomes exact in the limit of high spatial dimensions d or coordination number Z 
and may be employed to compute the phase diagram and thermodynamics of interacting lattice boson systems in the 
entire range of microscopic parameters. The B-DMFT requires a different scaling of the hopping amplitude with Z 
depending on whether the system is in the normal or the Bose-condensed phase. This additional difficulty compared 
to the fermionic case prevented the formulation of the B-DMFT in the past. As shown here it can be overcome by 
performing the scaling not in the Hamiltonian but in the action. The B-DMFT equations consist not only of a bosonic 
single-impurity problem in the presence of a self-consistency condition (the momentum integrated Dyson equation), 
but involve an additional coupling to the condensate wave function. 

We documented the comprehensive nature of the B-DMFT by explicitly reproducing results previously obtained in 
special limits of parameter space and by deriving other bosonic mean-field theories. For example, by calculating the 
local-self energy of the B-DMFT in the weak coupling regime, U/t <C 1, in perturbation theory and by including all 
normal and anomalous terms to first order in U/t one obtains the Hartree-Fock-Bogoliubov self-consistent mean-field 
approximation^ By neglecting the anomalous terms the standard Bogoliubov theory for lattice bosons is recovered. 
Inclusion of the second-order corrections to the local self-energy corresponds to the Beliaev-Popov approximation (with 
the additional assumption of a local self-energy). Furthermore, by neglecting all terms containing the hybridization 
function in the local action one obtains the mean- field theory developed in Ref. [TH and Refs. [l7lfl8l which corresponds 
to the exact solution of the bosonic Hubbard model |1| in the large dimension limit if only integer scaling is applied. 

In contrast to previous mean-field theories the B-DMFT constructed here treats normal and condensed bosons on 
equal footing. In particular, the B-DMFT takes into account effects due to finite hopping and dynamical broadening 
of the quantum levels. The inclusion of the hybridization function leads to important changes of the results of the 
static mean-field theory. For example, the Hubbard <5-peaks in the spectral functions now acquire a finite width, 
and the Mott transition occurs already when these bands start to overlap. We note that the B-DMFT is not merely 
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a perturbative improvement of the static mean-field theor y 13 i 17 i 18 i 20 with respect to the hybridization since in the 
B-DMFT the hybridization function is included to all orders. 

We applied the B-DMFT to solve the bosonic Falicov-Kimball model, i.e., a lattice model of itinerant and localized, 
interacting bosons. Due to the localized nature of the interacting bosons (providing a type of annealed disorder to the 
system since the localized particles are thermodynamically coupled to the itinerant bosons) the problem reduces to 
a set of algebraic equations. We find that local correlations enhance the transition temperature into the condensate 
and can thus be employed in the laboratory to increase Tbec- 

In general, the local single-site problem of a bosonic impurity coupled to two baths (the condensate and normal 
bosons) has to be solved numerically. The development of a reliable bosonic impurity solver— is a challenging task, 
which took several years in the case of the fermionic DMFT. This process can also involve a formulation of the proper 
bosonic impurity Hamiltonian corresponding to the B-DMFT action derived here. One of the main goals of this paper 
is to present the foundations of a novel, comprehensive mean-field theory for correlated bosons and thereby instigate 
further research by analytical and numerical means. 
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APPENDIX A: DERIVATION OF THE B-DMFT SELF-CONSISTENCY EQUATIONS 

Here we derive the B-DMFT equations for the generalized bosonic Hubbard Hamiltonian (JTJ) by applying the cavity 
method. 28 To this end the partition function 

Z = J D[b*Me X p(S[K,b v }) (Al) 

is calculated within the grand canonical ensemble, making use of the path integral approach over complex coherent 
states. 52 The action 

S[KM = / dT[J2bUr){d T -p,)b iu (T)+H(r)] (A2) 
Jo 



(A3) 



is split into a single-site term with i = 

S = / dT[b* Qv {T){d T - fi)b 0lJ {T) + - ^2 U^no^nau - S^)], 
a term representing the coupling between this site and the rest of the lattice (i ^ 0) 

AS = [ drJ2 (toiblhu + Cblbo u ) ee [ drAS(r), (A4) 
and a remaining part with site indices i,j =/= 

S<°> = [ P dr[J2 &*>X^-M)Mr) + ff (0) (r)L (A5) 

such that 

S[b* v ,b u ] = S + AS + S(°\ (A6) 

In the next step we expand the exponential function with respect to the action AS*, and perform the functional 
integral over all variables with site indices i 7^ 0. As a result we obtain a formally infinite series with all possible 
many-particle correlation functions, i.e., 

Z = J D[b* Ou ,b 0u \e- s °W»< b ^zM h - jyr{AS(T)) sm +^j\r 1 j\r 2 {AS(r 1 )AS{T 2 )} sl o ) +--A , (A7) 
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where (...)s(oj denotes the average taken with respect to (the action where the site i = is excluded) and 
is the corresponding partition function. In contrast to the fermionic case 2 ^ there remain anomalous correlation 
functions in the Bose-condensed phase such as (bi U (r)) S ( ) , {bi V {r)bj u {r')) g(o) , or (&* !y (r)&j ! y(r / )&fc 1 y(r")) S (o) , etc. The 
lowest first-order terms take the form 



/ dT(AS(r)) sm = [ P drJ2T! [«»<^ (t)) s(0) +^ b 0v (r){b* 1/ (T)) s(0) ] , (A8) 
Jo Jo v i 

where the prime on the summation symbol indicates that the lattice indices are different from 0, i.e. j ^ in (|A8[) . 
The second-order terms read 

-J d Tl J dr 2 (A^(r 1 )A5(r 2 )) 5( „ ) = -j dr x J dr^Y. 

u jk 

+to j tto( b Uri)bUr 2 )} s(0) bMb 0v (r 2 ) + ^^(^(n)6 fe ,(T 2 )) S (o,6S I/ (r 1 )6^(r 2 ) ] . (A9) 
Higher-order terms are obtained similarly. Defining the (connected) correlation functions for the condensate 

= (Mt)>sco) (A10) 
= (All) 

(|A8|) can be written as 

f dT(AS(T)) sm = f dr^zZ' [tojbUr^juir) + ^oM^-(t)] • ( Al2 ) 
Jo Jo v j 

Similarly, we define connected correlation functions for the one-particle excitations above the condensate as 

G]l (0) (ri - r 2 ) = -(T r MnK> 2 ))sc°>, (A13) 

GfJ 0) (n - r 2 ) = -(T r 6^(n)6fcv(75))s(o), (A14) 

GfJ 0) (n - r 2 ) = -{T r b ju (n)bk*(n))sm, (A15) 

G-L (0) (n - *s) = -(^6*„(n)6L(T 2 )) s ca). (A16) 
which permits us to express the second-order contribution, (IA9[) . as 

- / dri / dr 2 (A5(r 1 )A5(r 2 )) s(0) = -- / dn / dr 2 ^^ 
^ ! Jo Jo z - Jo Jo „ jk 

^o*feo G iL (0) ( T i - ^M'rOM'nO + t^okGfk^ (n - T 2 )6 „(ri)6o ;y (T 2 ) 

+ioj*fcoG J 1 L ( ° ) ( T i - t- 2 )^(ti)^(t 2 ) +toj*ok G ]L (0) ( T i " T 2)&L( T i)&o,.(' r 2) 
-iyfeo^v( T i)^L( T 2) & o^(-ri)^(T 2 ) - *5'o<ofe0i i /( T i)0fe 1 /(r 2 )&o ! ,(ri)&Si/( r 2) 
-*oi*fco^(n)0L(r a )65 v (ri)6o v (r a ) - WkMnWk^T^MJ&fo) ] . (A17) 

Here the first four terms are due to connected contributions and the last four terms due to disconnected contributions; 
higher terms can be written in a similar way. 

A non-trivial limit d — > oo is obtained by scaling the hopping amplitudes of i/-bosons as described in Sec. II. B. 
Namely, integer scaling is applied if t"j appears together with at least one anomalous average <fii V {T) = (bi V (r)) S (o) 
involving the BEC, while fractional scaling is employed otherwise. For example, in the first-order term, (|A8|) . the sum 
over j gives a contribution of the order 0(Z Roj ) so that the hopping amplitude <q ■ must be scaled with (i.e. divided 
by) a factor Z Roj because <pj V does not depend on the distance. On the other hand, in the first four terms of the 
second-order contribution to the partition function, (IA9|) . the hopping amplitudes must be scaled with Z R °^ 2 because 
the one-particle correlation functions are already proportional to l/Z Roj / 2 as discussed in the Sec. II. B. In the last 
four terms of the second order contribution the hopping amplitudes must be scaled with Z Roj . In the calculation of 
higher-order terms one has to distinguish the cases where all site indices are different from those where some, or all, 
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are the same. In analogy to the fermionic case, discussed in detail in Ref. l28|, we find that all connected higher-order 
terms vanish at least as 0(1/Z). Consequently, in the Z — > oo limit only connected contributions containing cj)j V 

or Gj%,£ \ or disconnected contributions made of products of connected contributions remain, provided the infinite 
series converges at least conditionally. Finally, we assume that the system is homogeneous, i.e., that <j>i V — <j) v is 
site independent. Applying the linked cluster theorem and collecting only connected contributions in the exponential 
function one obtains the local action 




«G-l ( Vi - r 2 )Mri)Mr 2 ) + % % k GfJ 0) (n - 7a)Mn)&S„fo) 
+% j %oG) 1 J 0) (.n - T 2 )bUri)b AT2) + *Vofc G }l ( °Vi - T 2 )b* 0v (n)b* 0v (T2) 

+^ U^no^(no v - Sftu), (A18) 

where the numerical factor k v = J2i^o^io/^ Ri " f° r ^ ~~ * 00 depends on the lattice structure. 

To simplify notations we introduce the Nambu formalism^ by defining a spinor boson operators hi V = (bw, b J „) and 
corresponding complex variables in the path-integral representation. Thereby anomalous averages for the condensate 

* iv (r) = (b w (r)) s(0 , (A19) 

and connected propagators for normal bosons 

Gj#(r - r') = -<T T b w (r)b] t/ (r')) S (o, (A20) 
can be written in a compact vector or matrix form. Introducing the hybridization matrix function 

A„(t - r') = -^%G|(r - r'), (A21) 

ij 

and employing the free ( "Weiss" ) mean- field propagator Q u one can express the B-DMFT local action in the form of 
(llOp . Here the site index i — is omitted for simplicity. 

Finally, the lattice self-consistency condition (fT2| needs to be derived. For this we apply the relation between the 
Green function G^(r — r') where the site i = is removed and the full lattice Green function, i.e. 

^*ijv = Giju — G^qi/Gqq^GojV) (A22) 

which holds for a general lattice. In the B-DMFT self-consistency equations ([ MT5)l for a homogeneous system only 
the site index i = enters which is therefore dropped. 

APPENDIX B: FREE BOSONS ON THE BETHE TREE WITH INFINITE COORDINATION NUMBER 

In this Appendix we employ the B-DMFT to study a single species (y = 1) of non-interacting bosons on the Bethe 
lattice with Z = ooi 54 i 55 Although this problem is exactly solvable by different methods^ it is instructive to see how 
the B-DMFT works in detail in this case. 



1. Green function method 

In order to obtain the Matsubara Green function 

G ij (r-r') = -{T T b i (r)b](r')) (Bl) 

for non-interacting bosons described by the Hamiltonian (TJJ we use the Bogoliubov transformation to separate the 
operator bi into a normal (non-condensate) part bi and the condensate wave function cf>{ as 
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hi = bi + fa 

b\=l\ + 4>l (B2) 

Assuming the system to be homogeneous, fa = <f), this yields 

G tJ (r - t') = -|</)| 2 + G l3 (r - t'), (B3) 

where Gij is the Green function of the normal bosons. In the non-interacting case considered here the anomalous 
Green function is absent. The density of particles is given by 

n = - lim ±J2 G ^ T - T ') = \^\ 2 - lim + ^E^( r - T ')' ( B4 ) 

z a 

where iVx, is the number of lattice sites. 

The diagonal Green function of normal bosons is given by 

G u (t -t') = ~J2 e-^ T -^G u {w n ), (B5) 

n 

where 

G u {u n ) = -Lj2~-^ r = ■ T ~~a7 v (B6) 

yVx, f-' «w„ + /i — A iw„ + /i - A(w„ ) 
and A are the exact energy eigenstates of the lattice Hamiltonian. The recursion relation 

GuM= . * , (B7) 

which is exact for the Bethe lattice ; 54 ' 55 allows one to express the hybridization function as 

AK)=t 2 G»K). (B8) 

Eq. (|B7|) determines Gjj as 

Gii(w n j = ^- . (B9j 

In particular, the equation for the particle density follows as 

n = \<j>\ 2 -±Y j / ( B1 °) 

Cauchy's theorem allows one to express the infinite sum as an integral over the spectral function multiplied by the 
Bose-Einstein distribution functio n 36 ' 52 ' 53 such that the density equation takes the form 



n = H 2 + -L [ ^ dco . U \ . (Bll) 

For temperatures T > Tbec the condensate vanishes, \ fa — 0, in which case the equation determines fi as a function 
of the density n. For T < Tbec the chemical potential is pinned at the value /i = — 2i and (|B11[) determines \<fi\ 2 , 
the density of the condensate. The condensation temperature Tbec itself is thus obtained for /i = — 2t and \(f>\ = 0. 
Expanding the Bose-Einstein function into a Taylor series and changing the integration variable into to — 2t cos 6 one 
obtains a transcendental equation for Tbec 

47rf2 °° -kx 

X K 
k=l 



where x = 2t/TBEC and I\ is a modified Bessel function 



57 
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2. B-DMFT 



We now show that the same results can be derived directly from the B-DMFT equations. The local action takes 
the explicit form 

rP rP rP 
S ioc = - / dr dT'b*{T)g-\T - T')b(T') +t dr (6*(r)0(r) + 4>*{T)b{r)) , (B13) 
Jo Jo Jq 

where the local Weiss Green function (an operator) is given by 

g-\T - t') = S(t - t') {-d T + ft) - A(r - r'). (B14) 

The hybridization function A(r) is determined self-consistently by eqs. (|llll2p . As in the fermionic case the relation 
A(cj„) = t 2 Q(uj n ) also holds for non- interacting bosons on the Bethe lattice in the limit Z — > oo. 

In the absence of interactions the Euler-Lagrange equation of motion for the classical field <p(r) is given by 



Q SS loc \b,b* 



6b* (r) 



= (8 T -IM)<f>(T)+ / dT'A{T-T')(f)(T')+t<f>{T). (B15) 







f,(r)=0(r) 

By Fourier transformation (|B15|) becomes a linear equation 

(i<jj n + [ji-i-A(w n ))<t>(w n ) = 0. (B16) 
Employing (|B9|B10p this equation takes the form 



iuj n + n-2t+ yj(iuj n + /i) 2 - (2t) 2 j <f>{w n ) = 0. (B17) 

In the static limit uj n = 0, corresponding to n = 0, this equation has the solution cj> = when fi < —2t, implying 
that the chemical potential lies outside the bosonic band, or the solution <p 7^ when /i = — 2t. In the latter case the 
actual value of \<p\ (which determines the BEC fraction) must be computed from the equation for the particle density 
(|Bllj) . In the dynamical case, w„ ^ 0, i.e. for u^O, (|B17|) only has the solution (j){u) n ) = because the expression 
in the bracket never vanishes. This shows that for non- interacting bosons the condensate order parameter is time 
independent. 

After Fourier transformation the local action takes the form 

Sioc = 6 » t^" + A 4 - A K)] b„ + t(b* n=0 (/) + b n=0 (j)*) , (B18) 

n 

where the numbers b n are the Fourier coefficients of 6(r) in Matsubara frequency space. The zero frequency component 
b n= o is shifted by the Bogoliubov transformation (|B2j) as 6 rl= o = b n =o + 4>- Because of (|B16|) the local action is seen 
to be quadratic in b n , and the functional integral yields the same equation for the particle density as (|B11I) . Thus we 
showed that the B-DMFT correctly reproduces all results for non-interacting bosons on the Bethe lattice both in the 
normal and in the BEC phase. 
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